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ABSTRACT 

Large-scale simulations of the Centaur population are carried out. The evolution of 
23 328 particles based on the orbits of 32 well-known Centaurs is followed for up to 
3Myr in the forward and backward direction under the influence of the 4 massive 
planets. The objects exhibit a rich variety of dynamical behaviour with half-lives 
ranging from 540 kyr (1996 AR20) to 32 Myr (2000 FZ53). The mean half-life of 
the entire sample of Centaurs is 2.7 Myr. The data are analyzed using a classification 
scheme based on the controlling planets at perihelion and aphelion, previously given in 
Horner et al (2003). Transfer probabilities are computed and show the main dynamical 
pathways of the Centaur population. The total number of Centaurs with diameters 
larger than 1 km is estimated as ^ 44300, assuming an inward flux of one new short- 
period comet every 200 yrs. The flux into the Centaur region from the Edgeworth- 
Kuiper belt is estimated to be 1 new object every 125 yrs. Finally, the flux from the 
Centaur region to Earth-crossing orbits is 1 new Earth-crosser every 880 yrs. 

Key vifords: minor planets, asteroids - planets and satellites: general - celestial 
mechanics, stellar dynamics - Kuiper belt 



1 INTRODUCTION 

The dynamical behaviour of the Centaurs is still poorly un- 
derstood. It is possible for a Centaur to work its way slowly 
inwards through the outer Solar system, leading to eventual 
capture by Jupiter and designation as a short-period comet. 
It is also possible for Centaurs to drift outwards to join the 
Edgeworth-Kuiper belt, to be ejected from the Solar system 
in an encounter with one of the massive, outer planets, or 
even to be captured by these planets into temporary satellite 
orbits. A small number may even impact upon the planets. 
Therefore, the Centaurs potentially hold the key to under- 
standing the mechanisms by which the short-period comet 
population is maintained, to explaining the distant, retro- 
grade satellites of the massive planets, and to allowing us 
a glimpse of objects newly introduced to the Solar system 
from the Edgeworth-Kuiper belt. 

As early as 1990, when the only known Centaur was Ch- 
iron, it was realised that such objects may lie on very unsta- 
ble orbits. Numerical integrations by Hahn & Bailey (1990) 
found that Chiron had a half-hfe for ejection of around 1 
Myr, but that the half-life to become a short-period comet 
for the object was around only 200 000 years, implying that 
Chiron could well have been a short-period comet in the 
past and could possibly become one in the future. This is of 



particular interest given the size of Chiron (d ~ 140 — 180 
km; Groussin, Lamy & Jorda 2004 and references therein) 
and other Centaurs, because objects that large entering the 
inner Solar system would be both spectacular and dan- 
gerous. In fact, the idea has been mooted that objects of 
such size arrive in the inner Solar system with some fre- 
quency, and then fragment, leading to swarms of debris 
which have the potential to encounter the Earth. The Kreutz 
sun-grazer family may represent one example of such hier- 
archical fragmentation, whereas other cases in which comets 
of more ordinary size have undergone catastrophic fragmen- 
tation include 3D/Biela, D/1994 (Shoemaker-Levy 9) and 
C/1994S4 (LINEAR). Whether such a decay mode repre- 
sents a generic process in determining the number of short- 
period comets can in principle be tested by examining the 
differences in the size distribution of Jupiter-family comets 
from those of their probable source objects, namely Cen- 
taurs, Edgeworth-Kuiper belt objects, long-period comets 
and so on (cf. Lowry, Fitzsimmons & CoUander-Brown 2003; 
Lamy et al. 2004). Clube & Napier (1984) have suggested 
that the Taurid meteoroid swarm may be the relic of the 
last large object to undergo such a decay. 

Despite the importance of the Centaurs, there has been 
little systematic study of the population using numerical 
simulations. Early calculations on Chiron (Oikawa & Ev- 
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Table 1. The classification scheme introduced by Horner et al. 
(2003). In the upper table, the first letter designates the planet 
controlling the perihelion, the second letter the planet controlling 

the aphelion or the region in which the aphelion lies, with the 
final two classes EK and T being beyond all the giant planets. 
(J = Jupiter, S = Saturn, U = Uranus, N = Neptune, EK = 
Edgeworth-Kuiper belt, T = trans-Neptunian belt). The lower 
table refers to cometary bodies (E = Encke, SP = short-period, 
I = intermediate and L = long-period). 



erhart 1979, Hahn & Bailey 1990) and Pholus (Asher & 

Steel 1993) identified the chaotic nature of these two ob- 
jects, tiiough only small numbers of clones and modest in- 
tegration times (< 1 Myr) were used. Work by Dones et 
al. (1996) looked at the behaviour of four Centaurs (Chi- 
ron, Pholus, Nessus and 1994 TA) and two Jupiter-family 
comets (29P/Schwassmann-Wachmann 1 and 39P/Otcrma). 
They found that the number of surviving objects decays ex- 
ponentially during the early part of the integrations, whilst 
the decay becomes flatter after a number of life-times. They 
also noted that Centaur half-lives inferred from numeri- 
cal integrations are smaller than those deduced from ap- 
proximations like Opik's (1976) theory and diffusion equa- 
tions (e.g., van Woerkom 1948). Lcvison & Duncan (1997) 
ran orbital integrations of 2200 test particles evolving from 
the Edgeworth-Kuiper belt to short-period comets, passing 
through the Centaur region in the process. The study of the 
integrations was mainly focused on the behaviour of the ob- 
jects both in the Edgeworth-Kuiper belt and the cometary 
region, rather than in the Centaur region itself. Hence, the 



dynamics of the Centaur population remains largely unex- 
plored. 

In an earlier paper, Horner et al. (2003) introduced a 
new method of classifying objects in the Solar system. This 
was based on the idea that the dynamical evolution is largely 
determined by the planets that control the perihelion and 
the aphelion of the object. This classification is particularly 
useful for the Centaurs, since it brealcs down the trans- 
Jovian region into 20 categories given in the upper panel 
of Table 1. Objects are labelled according to the control- 
ling giant planet, so, for example, a JS object has perihelion 
under the control of Jupiter and aphelion under the con- 
trol of Saturn. Objects with perihelion distance q > 33.5 au 
arc designated as either members of the Edgeworth-Kuiper 
belt (EK) or trans-Neptunian disk (T). In addition, objects 
with g < 4 au are designated as comets and are subdivided 
into Encke types (E), short-period (S), intermediate (I) and 
long-period (L), as summarised in the lower panel of Table 1. 

The aim of this paper is to provide results from a de- 
tailed set of simulations, exploiting the new classification 
scheme. The orbits of 32 of the best known Centaurs were 
used to create an ensemble of 23 328 clones. The clones were 
integrated in the presence of the four massive outer plan- 
ets in both forward and backward directions for a period of 
three million years, giving a vast data set with which to ex- 
amine the dynamics of their orbits. §2 describes the details 
of the numerical simulations, whilst §3 provides half-lives for 
individual Centaurs. §4 gives transition probabilities which 
allow the main dynamical pathways through this region of 
the Solax system to be identified. The simulations are used 
to estimate the total population of Centaurs in §5, together 
with the typical fluxes inwards from the Edgeworth-Kuiper 
belt. Finally, §6 considers possible correlations between the 
dynamics and the colours of the Centaurs. 



2 INTEGRATIONS 

In order to study the bulk statistics of Centaurs, 32 objects 
were selected from the list of Centaurs given on the Mi- 
nor Planet Center's website*. The objects were restricted 
to those with an observed arc of at least 30 days and an 
aphelion distance of less than 40 au. This ensures that only 
Centaurs with moderately well-determined orbits were in- 
cluded in our sample. The list of objects is given in Table 
2. Over time, as the Centaurs are observed over longer axes, 
the accuracy with which the orbits are known increases, and 
the orbits given on the Minor Planet website change accord- 
ingly. The orbits used in these integrations therefore repre- 
sent the best available information as of June 2002. Table 2 
also gives each Centaur's absolute magnitude H, which is 
defined as the apparent magnitude the object would have, 
if it were placed at both 1 au from the Earth and 1 au from 
the Sun and was observed at zero phase angle. This is calcu- 
lated in ignorance of any out-gassing that might occur. We 
can estimate the maximum and minimum diameter, assum- 
ing values of the albedo between 0.15 and 0.02. This gives 
a crude reckoning of the size, though photometric work is 
required to obtain any more detail. Of the objects studied 

* http://cfa-www.haj:vaj:d.edu/iau/lists/Centaurs.html 
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Table 2. The names of the objects simulated, arranged in order of increasing semi-major axis, together with their orbital elements as of 
June 2002. Here, a is the semi-major axis measured in astronomical units (au), e is the eccentricity, i is the inclination (in degrees), u) is 
the argument of perihelion (in degrees), f2 is the longitude of the ascending node of the orbit (in degrees) and H is the absolute visual 
magnitude. D^i^ and Dmax are the values for the diameter (in km) assuming albedos of 0.15 and 0.02 respectively, which are typical 
upper and lower limits from albedo measurements of comets and Centaurs performed to date. 'Class' is the classification of the object, 
using the scheme given in Horner et al. (2003). [The data are compiled from the Minor Planet Center.] 



© RAS, MNRAS 000, 1-15 



4 J. Horner, N. W. Evans & M.E. Bailey 



in these integrations, the one with the brightest absolute 
magnitude is 1995 SN55, with a value oi H = 6.0, which 
corresponds to a diameter of between 220 and 590 km. The 
object with the faintest absolute magnitude is 2000 GM137, 
with H = 14.3 giving a diameter of between 5 and 13 km. 
This is similar to the size determined for some cometary 
nuclei (e.g. Sanzovo et aJ. 2001; Lamy et al. 2004). Hence, 
the Centaurs come in a wide range of sizes, from very large 
(1995 SN55 and Chiron) to those comparable in size with 
normal comets (2000 GM137). 

The orbital elements of each object were used to cre- 
ate a swarm of 729 clones, distributed through a small cube 
of a-e-i (semimajor axis, eccentricity and inclination) space, 
centred on the original orbit. The clones of the objects were 
created by incrementally increasing (and decreasing) the 
semi-major axis of the object by 0.005 au, the eccentric- 
ity by 0.005, and the inclination by 0.01°. These increments 
arc sufficiently small that the clones can be considered as 
initially essentially indcntical to one another, yet they are 
large enough to ensure that the subsequent chaotic dynam- 
ical evolution following close planetary encounters rapidly 
disperse their orbits through phass space. So, nine values 
were used for each of those elements, with the central (fifth) 
value of the nine having the original orbital elements. This 
gives 729 clones of each of the 32 Centaurs, giving a grand 
total of 23 328 objects. 

The use of multiple clones of an individual object in the 
study of its behaviour over time is desirable for a number 
of reasons. First, the observations contain some uncertainty, 
which means that the orbit itself is not known beyond a 
certain degree of precision. This alone would bo enough to 
promote the use of a cluster of orbits with slightly differ- 
ent parameters. In addition, the chaotic nature of the orbits 
implies that an infinitesimally small change in the initial pa- 
rameters may lead to a major difference in the final outcome 
of the simulation. This moans that, beyond a certain time 
in the future, an object could be anywhere within the Solar 
system, or even beyond, as a result of a tiny change in the 
initial elements. These two facts taken together suggest that 
the best means to examine the future or past behaviour of 
an object is to integrate a large number of clones, and to 
examine the statistical properties of the dataset (Hahn & 
Bailey 1990, Dones et al. 1996). The number of clones used 
in such a simulation is chosen to maximise the size of the 
dataset available for analysis, without requiring an excessive 
amount of time for the simulations to run. The simulations 
described here took about three months to run on a desktop 
workstation. 

The clusters of Centaurs were then integrated for 3 Myr 
in both the forward and backward directions. The gravi- 
tational influence of the four Jovian planets (Jupiter, Sat- 
urn, Uranus and Neptune) was included in the integrations, 
which were all carried out using the hybrid integrator within 
the Mercury (Chambers 1999) software package. This is a 
symplectic integrator, which makes use of a turnover func- 
tion to switch to an accurate Bulirsch-Stocr algorithm for 
close encounters. The terrestrial planets (Mercury, Venus, 
Earth and Mars) were all omitted from the integrations, and 
their masses added to that of the Sun. The only slight detri- 
ment that this causes is loss of accuracy when objects are 
captured into orbits crossing those of the terrestrial planets. 
Even in this case, however, the effects of Jupiter (and the 



other giant planets, if the object's aphelion lies sufficiently 
far from the Sun) are generally much greater than those of 
the terrestrial planets. 

After some trials, a time step of 120 days was used for 
the integrations, since this was found to give a good com- 
promise between speed and accuracy. In order to determine 
the most efficient timestep, an object was placed on a typ- 
ical short period cometary orbit with perihelion near the 
Earth and aphelion near Jupiter. A number of clones were 
created, and the ensemble was integrated for 10'' yrs with 
time steps of 30, 60, 120, 240 and 360 days. The result- 
ing orbital elements were then compared, and it was found 
that the results for timesteps of 30, 60, 120 and 240 days 
gave consistent results, while 360 days was too long a time 
step. After a number of such trials, a time step of 120 days 
was used for the integrations, since this was found to give 
a good compromise between speed and accuracy, the mid- 
range value of 120 days being chosen so as to err on the side 
of accuracy where possible. 

An ejection distance of 1000 au from the Sun was 
used, following Lcvison & Duncan (1994). Any object which 
reached this distance was removed from the integration ^ . 
Also removed were those objects which impacted upon the 
surface of the Sun (g < 0.005 au), or on any one of the giant 
planets (the separation is less than the physical radius of 
the planet). On completion of the integration, datafiles were 
created for each clone which gave the values of the orbital 
elements at 100 year intervals for the clone's entire lifetime 
within the integration. It was on these files that the analysis 
was carried out. Each one of these files was approximately 
5 Mbytes in size (the exact size varied, since the file termi- 
nated with the ejection of the clone from the simulation). 
Hence, the 23 328 clones in total occupied « 120 Gbytes of 
disk space for their orbital elemental evolution alone, prior 
to any analysis - a daunting dataset by any standards! 



3 HALF-LIVES 

It is straightforward to calculate the value of the half-life 
for each Centaur. As illustrated in Figure 1, the number of 
clones remaining within the simulation decays in a roughly 
exponential manner as a function of time. The four objects 
whose decay is shown are 1996 AR20 (the object with the 
shortest lifetime of those studied), Pholus (an object with 
a moderately short lifetime), Nessus (a relatively long lived 
object) and 2000 FZ53 (the object with the longest lifetime). 
It has often been noted in long Solar system integrations 
that the number of clones remaining within a simulation 
decays exponentially with time (e.g., Dones et al. 1996, Hol- 
man 1997, Evans & Tabachnik 1999). The trait is less ob- 
vious with the longer lived objects. For example, Dones et 
al. (1996) found that the mirnber of surviving clones in their 
integrations decayed exponentially at early times, while at 
later times (generally greater than twice the half-life) the 
decay was slower. This is because those objects surviving 

t The long-period comets (the L class in Tabic 1) arc defined 
to have aphelion Q in excess of 1000 au. Our choice of ejection 
distance means that the statistical properties of the L class cannot 
be reliably computed from our simulations. 
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Figure 1. The number of surviving clones (lower left panel) and its natural logarithm (lower right panel) as a function of time for 1996 
AR20, Pholus, Nessus and 2000 FZ53. The upper panels show the same quantities but now for the entire ensemble of 23 328 clones. In 
all cases, the graphs are deduced from the simulations in the forward direction. 



the longest are the ones transferred to the most stable areas 
of the Solar system. 

The number of clones A'' which remain after time t is 
therefore given by 



N = Noe 



0.693 

Ti ' 



(1) 



where A'o is the initial number of clones, A is the decay 
constant, and Ti is the half-life for the object. To calculate 
the half-life, the simulation data were analysed with the help 
of least-square fitting routines from Press et al. (1992). The 
software provides the value of the function 



2,^^ [\n Ni - In No + \ti]^ 



1=1 



(2) 



where Ni is the number of clones remaining at time ti and 
the (Ti are the individual standard deviations on the data 
points. As the ai are unknown, we proceed by first assigning 
uniform errors, fitting for the model parameters by minimiz- 
ing the ™d then rescaling the errors using eqn. [15.1.6] 
of Press et al. (1992). Of course, this well-known procedure 
precludes an independent estimate of the goodness of fit. 

The overall dataset of 23 328 objects has an ensemble 
half-life of 2.76 Myr in the forward direction and 2.73 Myr 
in the backward direction. This gives us an estimate of the 
mean lifetime of a typical Centaur. Note that this lifetime 
adds weight to the argument that the population of the Cen- 
taur region is in a steady state (a reservoir of objects con- 
stantly being drained by Jupiter, and refilled from a long- 
lived source like the Edgeworth-Kuiper belt). 

Dones et al. (1996) calculated the half- lives of Chiron, 
Pholus, Nessus and 1994 TA using ~ 100 clones. Their re- 
sults for Chiron and Nessus are in excellent agreement with 
ours, but they found Ti — 2.1 and 2.4 Myr for Pholus and 
1994 TA - somewhat larger than our results. The most likely 
cause of the discrepancy is in the different algorithms used 
to populate the clones. Dones et al. carried out a compari- 
son with approximations like Opik's (1976) theory and dif- 
fusion equations (e.g., van Woerkom 1948) and concluded 
that both methods significantly overestimate the lifetimes 
by factors between 2 and 5. It seems that numerical simula- 



tions with large numbers of clones are the only reasonably 
reliable method for half- life estimation. 

The half-lives of individual Centaurs are given in Ta- 
ble 3. The value of the Poisson uncertainty a is calculated 
as a — Ti/2/\/ No, where A^o is the initial number of clones. 
This uncertainty is added in quadrature to the uncertainty 
in the fitted half-life, as judged from the surface in the 
space of model parameters and as returned by our fitting 
software. Note that our ignorance of individual error bars 
on our simulation datapoints may lead to an underestimate 
of the latter quantity - it is generally smaller than a by 
an order of magnitude - as our algorithm is tantamount to 
assuming a good fit to the exponetial decay law. The ob- 
ject with the shortest half-life is 1996 AR20, a JN object, 
which has a half-life of approximately 540 kyr in the for- 
ward integration and 590 kyr in the backward integration. 
The object with the longest half-life is 2000 FZ53, a UE 
object, with half-lives of 26.8 Myr (forwards) and 32.3 Myr 
(backwards) . 

On comparison with the orbital elements (Table 2), a 
correlation can be seen between the position of a Centaur's 
orbit within the Solar system and its half-life. The farther 
from the Edgeworth-Kuiper Belt, the shorter the half-life. 
This is not unexpected - Jupiter is significantly more mas- 
sive than Saturn, and Saturn in turn is more massive than ei- 
ther Uranus or Neptune. So, the farther from the Edgeworth- 
Kuiper Belt, the more massive are the planets with which 
the Centaur interacts and the more frequently do such en- 
counters occur. This effect is evident when the orbital el- 
ements of the objects are plotted against the logarithm of 
the half-life, as in Figure 2. There is only a rough correlation 
with semi-major axis, but the data indicate a lower bound 
to the half-life as a function of perihelion q and an upper 
bound as a function of aphelion Q. Specifically, we find that 



0.392 exp(0.135g) < T1/2 < 0.064 exp(0.275 Q). 



(3) 



where q and Q are in an and T1/2 is in Myr. This holds for 
all the Centaurs in our sample, but it is conceivably possi- 
ble that low eccentricity orbits between the planets are ex- 
tremely long-lived (e.g., Holman 1997; Evans & Tabachnik 
1999). Additionally, in the plot of perihelion distance versus 
half-life, there seem to be three rough groupings of objects. 
The first are those in a band along the line from q — 6, 
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Table 3. Half-Lives Ti (in Mvr) of the simulated objects, together witii an error estimate a (in Myr). The column labelled C gives the 

2 

class of the Centaur, while the column labelled D gives the direction of integration (F for forward, B for backward). 



logTi = 5.75 to q = 26, logTi = 7.2, which accounts for 
the bulk of the objects. A second group comprises six ob- 
jects, which lie roughly on a parallel track at values of log Ti 
greater by 0.7. Finally, 2000 FZ53 sits alone far above either 
of these groups. The objects which are in the first group cur- 
rently lie away from the positions of any major resonances 
and so tend to have short lifetimes. The objects in the sec- 
ond group tend to lie nearer to stable mean motion or sec- 
ular resonances. In fact, the orbit of 2002 GBIO lies within 
0.009 au of the 3:4 resonance of Uranus, while the orbit of 
Chariklo lies within 0.09 au of the 4:3 resonance of Uranus. 
Clones of 2000 FZ53 quite frequently display resonant be- 
haviour during the course of the simulations, although it 
does not currently lie near any major mean motion reso- 



nances. A possible cause of 2000 FZ53's exceptionally long 
half-life is its abnormally large inclination - higher than that 
of any other Centaur studied by over 10 degrees. Any cor- 
relations of half-life with inclination and ecccentricity arc 
less clear-cut than those with position. However, there is a 
lack of long-lived objects at large e. This is a consequence 
of the fact that Centaurs with large e must cross the orbits 
of several of the outer planets, and so inevitably are more 
unstable than bodies whose close approaches are restricted 
to just one or two planets. 

As the equations of motion are time- reversible, it might 
naively be expected that the forward and backwards half- 
lives should be the same. In fact, it is often found that inward 
evolution of the orbits of minor bodies is more likely than 
outward evolution (e.g., Oikawa & Everhart 1979). This is 
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Figure 2. The relationship of the logarithm of the half-lives of the Centaurs with their semi-major axis, perihelion and aphelion distances, 
and their inclination. Points in red show the results of the forward direction, and those in green show the results of the backward direction 
integrations. 



because the Edgeworth-Kuiper belt provides a source while 
Jupiter provides a sink, so that clones ejected by Jupiter 
do not have an opportunity to return. Of all the objects, 
19 have shorter half-lives in the forward direction, and 13 
have shorter half-lives in the backward direction. The least 
discrepancy occurs for 2001 BL41, with the two half-lives 
agreeing to within 0.25%. The greatest discrepancy occurs 
for 2002 GZ32, where the backward half-life is 46% longer 
than the forward half- life. These differences can be visualised 
as the effect of the first encounters the clones have with the 
major planets. If an encounter with one of the planets occurs 
early enough in the simulation, the entire ensemble of clones 
will be perturbed, moving the objects onto slightly differ- 
ent orbits, with slightly different half-lives. For the entire 
dataset, the forward and backward half-lives only diverge 
by a matter of 30000 years out of 2.76 Myr - a discrepancy 
of only just over 1%. 

From the dataset, it is also possible to calculate the 
half-lives of the starting class of the objects. The results of 
this calculation are given in Table 4. The number of clones 
in any particular class is not necessarily an exact multi- 
ple of 729. This is because a number of the objects have 
outlying clones which actually fall into a different class at 
the start of the integration as compared to the seed. The 
results of this analysis again show the dependence of half- 
life on perihelion position - objects in the Jupiter classes 
have half-lives shorter than those in the Saturn classes, and 
these in turn are more short-lived than the objects in the 



Uranus classes. There are also hints in the table that more 
eccentric objects under the control of any particular planet 
may be more long-lived than their less eccentric counter- 
parts (compare, for example, the half-lives of U, UN and 
UE objects). Orbital periods of the more eccentric objects 
are greater than those on near-circular orbits with similar 
perihelion distances, and hence encounters with the giant 
planets happen less frequently. 

Over the course of the integrations, the clones of each 
object are repeatedly transferred between classification bins. 
This allows us to evaluate the amount of time that is spent 
in each of these classes over the simulation, together with 
the number of times the object is transferred into that class. 
From this, we can calculate the mean time that an object 
spends in any particular class before being transferred into 
another. The results of such calculations are presented in 
Table 5. In this table, the value of the mean lifetime for 
the L or long-period comet class has been ignored, since 
objects which enter this classification are then removed from 
the simulation as they pass 1000 au from the Sun. This 
means that the value of mean lifetime for objects of class L is 
unrealistically small. It is also noteworthy that the EK and T 
classes have particularly short mean lifetimes. A stable orbit 
in these regions requires decoupling from Neptune, and there 
are no non-gravitational forces within the integrations which 
could allow this to happen. Hence, the very small number 
of objects which attain these two classes only do so at the 



© RAS, MNRAS 000, 1-15 



8 J. Horner, N. W. Evans & M.E. Bailey 



c 


Direction 


N 


Ti 

2 


a 


JS 


Forward 
Backward 


729 
729 


0.67 
0.65 


0.03 
0.02 


JU 


Forward 
Backward 


729 
729 


0.61 
0.64 


0.02 
0.02 


JN 


Forward 
Backward 


891 
891 


0.59 
0.63 


0.02 
0.02 


S 


Forward 
Backward 


729 
729 


0.72 
0.68 


0.03 
0.03 


SU 


Forward 
Backward 


3519 
3519 


1.03 
1.02 


0.02 

0.02 


SN 


Forward 
Backward 


2205 
2205 


1.11 
1.22 


0.02 
0.03 


SE 


Forward 
Backward 


1917 
1917 


1.75 
1.77 


0.04 

0.04 


U 


Forward 
Backward 


3780 
3780 


3.72 
3.36 


0.06 
0.06 


UN 


Forward 
Backward 


5094 
5094 


5.84 
6.19 


0.08 
0.09 


UE 


Forward 
Backward 


3006 
3006 


12.5 
14.2 


0.23 
0.26 


N 


Forward 
Backward 


729 
729 


12.5 
13.5 


0.46 
0.50 



Table 4. The half-lives of the individual classification bins across 
which the objects fall at the start of the integrations. N is the 
number of clones in that particular class at the start of the inte- 
grations, Ti is the half-life (in Myr), and <T is again the uncer- 

2 

tainty on the half-life (in Myr). 

extremes of a series of perturbations and are immediately 
perturbed back into classes under the control of Neptune. 



4 TRANSFER PROBABILITIES 

It is also straightforward to calculate the probability of an 
object being transferred from one cleiss to another. This can 
be visualised by constructing a 24 x 24 grid with the initial 
class on the vertical axis and the final class on the horizon- 
tal (these classes are just initial and final with respect to a 
single transfer, not for the entire integration). This is done 
by recording every transfer which occurs within the integra- 
tions, and hence calculating the fraction of objects which, for 
example, are transferred from class J to class JS. The results 
are shown in Table 6. The numbers have been normalised 
so that the sum along any row is unity. For any class, the 
probabilities give the relative likelihood of leaving from that 
class to the target classes given on the horizontal axis. As 
an example, let us take a typical result from one of these 
tables, namely that the value of the probability of transfer 
from class J to class JS is ~ 0.49 (Table 6). This means that, 
for an object in the J class, there is a ~ 49% chance of the 
object being transferred directly into the JS class the next 
time the classification changes. We can see that for such an 
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Table 5. The mean time (in yrs) spent in each classification 
bin, before a Centaur clone is transferred to another bin. The 
shortness of this mean time is understandable, as objects close 
to the boundary are often tranferred to and fro. Note that EC 
stands for Earth-crossing objects and MC for Mars-crossing. 

object, the two most likely transfers are to the JS class or 
to a short-period (SP) comet, and between them, these two 
possibilities make up the great bulk of transfers for all J 
objects. 

Table 6 shows a number of interesting features. When- 
ever an object is controlled by two planets (one at perihelion 
and one at aphelion) , the classes to which the object is most 
likely to move correspond to transitions at the perihelion 
and aphelion of the planet. For instance, an SU object is con- 
trolled at perihelion by Saturn and at aphelion by Uranus. It 
is most likely to be transferred to one of the classes JU or U 
by an encounter at aphelion. These cases corresponds to an 
encounter with Uranus either increasing the eccentricity of 
the orbit, and hence pushing the perihelion down to Jupiter's 
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control, or decreasing the eccentricity, pulling the perihelion 
away from Saturn's control. For encounters at perihelion, 
the most likely classes are S or SN corresponding either to a 
circularisation of the orbit at Saturn, or to a pumping of the 
eccentricity of the orbit, as the aphelion moves from Urarms' 
to Neptune's control. These most popular transfers can be 
traced diagonally down the tables, around the empty diag- 
onal corresponding to the same initial as final class. These 
four parallel lines of high probabilities give the appearance 
of two sets of "tram lines" running through the tables. Af- 
ter these possibilities, other transfers are also viable, albeit 
with lower probabilities - for example, an SU object can suf- 
fer a perihelion-aphelion interchange at Saturn, moving to 
the JS class. However, the fact that the four classes most 
likely to be reached in a transfer lie along the "tram lines" 
vindicates the classification scheme, which is based on the 
idea of transfers by interaction primarily at perihelion and 
aphelion. 

Also of interest is an effect which can be seen on com- 
paring Table 6 with the equivalent results for each individual 
Centaur (given in Appendix A of Horner 2003). Prom any 
class, the probability of transfer to another class is roughly 
constant, regardless of the direction of integration or the 
object in question. The main discrepancies lie in very low 
probability transfers, where the uncertainty is large because 
of the small numbers involved. This means that for a newly 
discovered object, it is possible to give the probabilities of 
its transfer to any new class, as long as the initial class can 
be computed. It also permits insight into the main dynami- 
cal pathways followed by a Centaur. For example, using the 
values in Table 6 and assuming an initial population of 1000 
short-period comets or SP objects, it can be seen that 27 of 
these objects become E types, all of which would return to 
the SP bin, 41 enter the I class, three of which would on av- 
erage then return, 304 travel to the J class, 147 making the 
return trip. For the other bins, we find 84 returning objects 
from JS, 8 from the JU class and one from ,JN. Therefore, 
270 of 1000 objects that leave the short-period class return 
immediately the next time that their classification changes. 

Given that the total time spent in each class and the 
number of times that class becomes occupied are calculable 
from the simulation data, it is also possible to compute the 
probability per unit time of a transfer. Wc already know 
the probabilities that an object in a particular class will be 
transferred to any other. Dividing the mean time spent in 
any class by this probability, we obtain the moan transfer 
time from one class to any target class. The inverse of this 
is the probability per unit time of the transfer. The values 
obtained in this type of analysis are given in Horner (2003) 
for the individual Centaurs, while the probabilities for the 
entire dataset are given in Table 7. The results are given as 
probabilities per Myrs, so that a value of 0.1 in a particular 
box means that an object making the relevant transfer would 
have a mean transfer time of about 10 Myrs. This means that 
the population in any class N evolves according to 

i 3 

where Pj and A^j are the probabilities per unit time and the 
populations in the bins along the row, while the Pi and Ni 
are the corresponding quantities along the column. In other 
words, the ingress to a particular class is governed by the 



numbers in the column, and the egress by the row. Math- 
ematically speaking, this gives us coupled sets of exactly 
solvable linear first-order difiierential equations that govern 
the evolution of Centaur clones. We will return to this in a 
later publication, but these ideas are already prefigured in 
Bailey et al. (1992). 



5 CENTAUR FLUXES AND POPULATION 

During the simulations, we also record the numbers of clones 
that become Earth-crossing objects. Mars-crossing objects 
and short-period comets. This gives us a moans to reckon the 
total population of Centaurs. Fernandez (1985) suggested a 
flux of 10^'^ new short-period comets per year, with a mean 
lifetime of ~ 6 kyr. More recently, Levison & Duncan (1994) 
find that the mean fade time for a short-period comet is ~ 12 
kyr. Under the assumption that the current population of 
short-period comets is in steady state, the work of Levison 

6 Duncan implies a flux of ~ 0.5 x 10~^ new short-period 
comets per year. This is equivalent to one new short-period 
comet being captured, on average, every 200 years. If we 
assume that the entirety of this flux comes from the Centaur 
region, then this allows us to estimate the total population 
of the Centaur region using the simulations. 

From the simulation data, a total of 7 900 out of 23 328 
clones 34% of the initial population) become short-period 
comets at some point during the forward integrations, and 
8 068 (again ~ 34%) become short-period comets during the 
backward integrations. This is a flux of one new short-period 
comet every 380 years. If we assume all short-period comets 
are captured from the Centaur region, an estimate of the 
total number of Centaurs (with perihelion distance g > 4 
au and aphelion distance Q < 60 au) is 44 300. This rep- 
resents the population of objects bright enough to be vis- 
ible as short-period comets, were they to be captured into 
such an orbit. An effective nuclear diameter d greater than 
1 km seems a reasonable limit to place for this value, though 
there are an increasing number of comets with d ~ 0.5 km 
(e.g., Lamy et al. 2004). This calculation also only takes 
into account objects becoming short-period comets for the 
first time. Given that the mean fade time is ~ 12 kyr, it 
is reasonable to assume that objects captured for the first 
time axe significantly brighter as short-period comets than 
those which have experienced a number of prolonged stays 
in the region. Therefore, the objects, which in the simula- 
tions display a number of prolonged periods as short-period 
comets, would actually exhaust all their volatiles early on 
and should not contribute to the new short-period comet 
flux in later passages through the inner Solar system. Note 
that the usage of the flux of short-period comets to normalise 
the source populations has also been exploited recently by 
Emel'yanenko, Asher & Bailey (2004) to estimate the to- 
tal population of trans-Neptunian objects on highly eccen- 
tric orbits, which reside at still greater heliocentric distances 
than the Centaurs. 

Given our estimate of the total population of Centaurs 
and the knowledge that the half-life is ~ 2.75 Myr, we can 
estimate the influx of now Centaurs from the Edgeworth- 
Kuiper belt. Neglecting those few objects that could be 
captured onto Centaur-like orbits from high eccentricity or- 
bits from the Oort cloud, we can see that ~ 22150 objects 
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0.253 






0.004 


0.360 


0.006 
















JT 






0.063 


0.009 








0.001 


0.055 








0.001 


0.002 


0.868 
















S 




0.001 








0.349 


0.006 










0.641 


0.001 




















SU 












0.006 


0.094 


0.001 






0.136 




0.449 


0.001 




0.311 


0.001 












SN 














0.002 


0.041 


0.001 






0.269 




0.312 




0.003 


0.366 


0.005 










SE 
















0.001 


0.076 


0.001 




0.001 


0.292 




0.152 




0.001 


0.475 


0.001 








ST 








0.010 










0.001 


0.563 








0.111 








0.001 


0.314 








U 
























0.330 


0.005 








0.664 












UN 
























0.001 


0.212 


0.001 




0.218 




0.449 




0.117 


0.002 




UE 


























0.003 


0.304 


0.001 




0.470 




0.061 


0.001 


0.159 




UT 




























0.001 


0.701 






0.159 






0.002 


0.137 


N 


































0..501 


0.005 






0.494 




NE 


































0.003 


0.492 


0.002 


0.374 




0.125 


NT 




































0.001 


0.566 




0.432 





Table 6. Transfer probabilities for the entire simulation in the forward (upper table) and backward (lower table) direction. Notice that 
the probabilities are nearly the same regardless of the direction of integration. An empty entry means that the transfer probability is 
< 10~^. The leading diagonal is empty by definition. 



must be replaced every 2.75 Myr. This is equivalent to 1 ob- 
ject transferred to the Centaur region from the Edgcwortli- 
Kuiper belt every ^ 125 yrs. This calculation also ignores 
the small flux of objects from Centaur-type orbits to the 
Edgeworth-Kuiper belt. Our simulations allow neither for 



non-gravitational effects, such as collisions, nor for the grav- 
itational perturbations between Edgeworth-Kuiper belt ob- 
jects. So, it is impossible to determine the flow of objects 
from orbits which encounter Neptune to those which are sta- 
ble beyond. To drive the inward flux of Centaurs, the efiiects 
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Pre 


Post 














































E 


SP 


I 


I 


.1 


.IS 


JU 


JN 


JE 


JT 


s 


SU 


SN 


SE 


ST 


U 


UN 


UE 


UT 


N 


NE 


NT 


E 




247.24 


0.15 










0.31 








0.23 


0.23 


0.08 




0.39 


0.31 












SP 


16.19 




24.33 


0.10 


181.12 


251.13 


86.67 


27.69 


5.24 


1.29 


0.29 


0.57 


0.32 


0.20 


0.02 


0.22 


0.17 


0.07 










I 


0.06 


122.96 




9..54 




0.08 




0.74 


144.25 


14.54.53 




0.63 


0.74 


0.74 


6.65 


0.85 


0.25 


0.11 










L 






1.100 






4.89 




21.70 




326.64 




51.42 


45.14 


26.68 


627.59 


46.19 


16.45 


13,56 


10.45 






0.11 


J 




1661.87 




0.09 




1711.44 


53.36 


3.95 


1.88 


0..50 


0..iO 


0.23 


0.32 


0.18 




0.05 


0.05 












.IS 




220.88 




0.02 


160.84 




521.23 


3.96 


1..35 


0.44 


193.99 


11.45 


0..34 


0.19 


0.03 


0.04 


0.02 


0.03 










JU 




42,.i7 




0.01 


3.02 


290.06 




260.57 


5.48 


1.10 


1.84 


1.50.93 


4.40 


0.28 


0.10 


0.05 


0,02 


0.03 










.IN 




19.51 


0.02 


0.03 


0.21 


2.47 


375.78 




323.20 


4.00 




2.98 


170.75 


5.76 


0.14 


0.06 


0.02 


0.02 


0.01 








JE 




2.62 


14.83 


0.10 


0.05 


0.50 


4.75 


231.08 




170.39 




0.10 


2.76 


248.55 


3.88 


0.09 


0.06 


0.03 










JT 




0.08 


60.54 


10.71 


0.01 


0.10 


0.24 


0.92 


61.89 






0.50 


0.63 


1.77 


1063.87 


0.38 


0.18 


0.10 










S 




0.87 






0.16 


260.30 


4.39 


0.01 


0.01 






450.08 


0.79 


0.20 


0.04 


0.01 


0.01 












SU 




0.11 








1.62 


26.65 


0.34 






36.99 




118.60 


0.35 


0.07 


87.64 


0.46 


0.01 










SN 




0.02 








0.06 


0.6.) 


16.48 


0..33 




0.08 


00.38 




118.33 


0.15 


1.23 


130..33 


2.22 










SE 




0.01 








0.03 


0.03 


0.42 


25.23 


0.37 


0.01 


0.24 


94.71 




49.43 


0.01 


0.48 


156.10 


0.17 








ST 






0.0.) 


3.07 




0.01 


0.01 


0.01 


0.28 


204.61 




0.09 


0.11 


36.76 




0.03 


0.04 


0.28 


92.48 








U 
























50.49 


0.85 


0.01 






98.95 


0.02 










UN 












0.01 












0.15 


41.35 


0.19 




45.80 




96.93 




21.C'S 


0.,32 




UE 












0.01 












0.04 


0.75 


64.34 


0.14 


0.01 


102.28 




12.46 


0.22 


,32,16 


0.04 


UT 








0.18 




0.04 




0.01 








0.12 


0.05 


0.56 


430.55 


0.01 


0.01 


106.75 






1.63 


99.46 


N 












0.03 












0.02 


0.01 


0.02 






187.82 


1.94 






195.65 


0.01 


NE 












0.05 












0.03 


0.01 


0.04 






1.08 


107.07 


0.64 


77.16 




29.37 


NT 












0.11 












0.05 


0.04 


0.07 






0.01 


0.38 


138.48 


0.01 


104.97 





Pre 


Posi 














































E 


SP 


I 


I, 


J 


JS 


jr 


JN 


JE 


JT 


s 


SU 


SN 


SE 


ST 


11 


UN 


TIE 


UT 


N 


NE 


NT 


E 




284.75 


0.54 








0.18 




0.09 




0.09 


0.54 


0.18 


0.27 




0.09 


0.27 


0.09 










SP 


16.43 




23.SS 


0.11 


189..36 


259.33 


90.64 


27.20 


3.,30 


1.,33 


0.26 


0.,3S 


0.,39 


0.22 


0.03 


0.04 


0.14 


0.11 


0.01 


0.01 






I 


0.05 


94.91 




8.36 


0.02 




0.55 


0.21 


109.06 


1372,92 


0.23 


0,46 


0,80 


0.48 


6.08 


0.16 


0.27 


0.30 


0.02 


0.02 






L 




4.,56 


93.76 




4.93 


9.11 


140.07 


39.48 


36.06 


210C,71 


1,34,76 


213.72 


311.65 


214.09 


4424.25 


113.88 


171.58 


107.05 


61 .,50 


7.97 


6.83 


0.76 


J 




1641.73 








1697.47 


55.66 


2.78 


1.12 


1,21 


0,81 


0.18 


0.18 


0.13 




0.04 


0.04 












JS 




224.80 




0.01 


166.55 




526.67 


3.90 


1.33 


0,56 


182,86 


11.37 


0.25 


0.14 


0.04 


0.02 


0.04 


0.04 










JU 




43.74 


0.01 


0.04 


2.72 


293.69 




264.54 


5.65 


1.21 


1.89 


146.54 


4.09 


0.23 


0.10 


0.02 


0.03 


0.03 




0.01 






JN 




19.99 


0.02 


0.06 


0.19 


2..)0 


384.13 




331.39 


3.86 


0.03 


2.56 


165.74 


5.66 


0.14 


0.02 


0.06 


0.03 










JE 






13.92 


0.11 


0.05 


0.49 


4.84 


235.46 




172,20 


0.04 


0.06 


2.61 


243.12 


3.96 


0.04 


0.06 


0.04 










JT 




0.14 


72.87 


10.80 


0.01 


0.08 


0.50 


0.76 


64.19 




0.21 


0.35 


0.65 


1.82 


1009.83 


0.18 


0.24 


0l14 










S 




0.84 






0.21 


258.23 


4.56 


0.02 




0,01 




474.09 


0.84 


0.22 


0.06 


0.01 


0.01 


0.02 










SU 




0.09 






0.01 


1.58 


25.99 


0.28 






37.42 




123.64 


0.36 


0.07 


85.52 


0.40 


0.03 




0.01 






SN 




0.03 








0.06 


0.61 


15.18 


0.29 




0.08 


99.24 




115.01 


0.14 


1.08 


135.05 


1.96 










SE 




0,01 








0.02 


0.04 


0.41 


25.10 


0.36 


0.02 


0.23 


96.90 




50.43 


0.01 


0.49 


157.74 


0.23 








ST 




0,01 


0.05 


,3.11 






0.02 




0.29 


182.06 


0.02 


0.06 


0.13 


35.84 




0.04 


0.03 


0.34 


101.61 








U 














0.01 










53.71 


0.85 








108.00 


0.05 










FN 














0.01 








0.01 


0.12 


42.19 


0.17 




43.21 




89.22 




23.31 


0.32 




i:e 














0.02 








0.01 


0.01 


0.70 


62.31 


0.19 


0.01 


96.,50 




12.54 


0.23 


32.71 


0.03 


UT 








0.15 






0.07 








0.04 


0.03 


0.13 


0.75 


474.46 




0.08 


107.39 






1.11 


92.84 


N 














0.01 








0.01 




0.02 








200.31 


1.84 






197.29 


0.01 


NE 




0,01 










0.02 








0.03 


0.01 


0.05 




0.01 




1.15 


113.18 


0.45 


86.08 




28.74 


NT 




0.02 






0.01 




0.04 








0.07 


0.01 


0.09 


0.02 


0.01 




0.03 


0.24 


127.54 


0.01 


97.38 





Table 7. Transfer probabilities per Myr for the entire simulation in the forward (upper table) and backward (lower table) direction. An 
empty entry means that the transfer probability per Myr is < 10~^. The leading diagonal is empty by definition. 



of collisions and of perturbations between the Edgeworth- 
Kuiper belt objects must be considered. Durda & Stern 
(2000) suggest that collisions of objects greater than 4 m 
in diameter onto comet-sized bodies within the Edgeworth- 
Kuiper belt occur every few days. In addition, they reckon 



that the time scale for the disruption of 1 km objects is ~ 1 
Gyr. These two facts in concert imply a high rate of collision 
within the Edgoworth-Kuiper belt, sufficiently high that an 
inward flux of 1 new Centaur every ~ 125 yrs seems reason- 
able. 
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From the forward integrations, we can deduce that 1799 
clones (7.7% of the sample) became Earth-crossing and 3799 
clones (16%) became Mars-crossing. Very similar numbers 
are yielded by the backward integrations. Therefore, we ex- 
pect typically one Centaur to become Earth-crossing for 
the first time approximately every ~ 880 yrs, and one new 
Mars-crosser every ~ 420 yrs. Most of the known popula- 
tion of Near-Earth objects (NEOs) is asteroidal in nature, 
as the Main Belt provides the great majority of NEOs. Mor- 
bidelli et al. (2002) state that only 6% of the NEOs arc ulti- 
mately of Edgeworth-Kuiper Belt origin. NEOs originating 
in the Main Belt are expected to survive for far longer times 
within the inner Solar system than the comotary bodies, 
due to the fact that they lie on orbits significantly decou- 
pled from direct perturbations by Jupiter. An approximate 
flux of one new Earth-crosser of Centaur (and hence, orig- 
inally Edgeworth-Kuiper belt) origin per millennium seems 
to be in reasonable agreement with such work. 

Though the number of NEOs of Edgeworth-Kuiper belt 
origin is small, they are exceptionally important in judging 
potential hazards. The size distribution of NEOs of aster- 
oidal origin is heavily weighted towards small particles, con- 
sistent with the idea that they arc collision fragments. The 
largest NEO is (1036) Ganymed with an absolute magnitude 
H = 9.45 and a diameter d w 50 km. Very few NEOs are 
larger than 10 km across. However, for the Centaurs, the 
upper end of the size distribution is well-populated, with 16 
objects of the 32 listed in Table 2 intrinsically brighter than 
Ganymed. The passage of a large Centaur like Chiron or 
Pholus into the inner Solar system would provide a very sig- 
nificant environmental disturbance (Hahn & Bailey 1990), 
as its fragmentation and possible decay could overwhelm 
the local space environment with debris and dust. 

Another interesting set of statistics is the total number 
of times a class becomes occupied during the integrations. 
In the forward direction, the Earth-crosscrs become occu- 
pied ~ 23 300 times over three million years. However, we 
already know that only 1799 of the clones enter this class, 
so it is obvious that the bulk of objects which become Earth 
crossing at least once will in fact enter, leave and re-enter 
the class a number of times (in this case, ~ 13 times). Ta- 
ble 8 summarises the number of visits to the short-period 
(SP), Earth-crossing (EC), Mars-crossing (MC) and Encke- 
type (E) classes that a typical object entering these classes 
would have. 

The flux of objects into the Enckc class of comets can 
also be calculated. This is subject to greater errors - be- 
cause of the much lower numbers of objects that become 
Encke-type comets than become Earth-crosscrs and because 
a time step of 120 days is insufficient to track the orbits 
of such objects with reasonable accuracy and because the 
gravitational influence of the neglected terrestrial planets is 
now significant. For objects in the forward integrations, 303 
clones became Enckc types, implying a capture rate of one 
every ~ 5200 yrs. The backward integration yields the same 
capture rate. 

Typically, an object will become an Earth-crosser or 
Mars-crosser while classified as a short-period comet, having 
been transferred to the inner Solar system by an encounter 
with Jupiter. Although some objects can be perturbed onto 
Earth and Mars crossing orbits by the other outer planets 
(for example, the comets P/1997 T3 Lagerkvist-Carsenty 



Class 


D 


N 


EC 


F 


12.9 






1.S.7 


MC 




1 7 




B 


14.6 


SP 


F 


28.0 




B 


28.4 


E 


F 


16.5 




B 


18.3 



Table 8. Table showing tiic mean number of times an object 
which enters a conictary class at least once will go on to enter and 
re-enter that class through its lifetime. D gives the direction of the 
integration (F = Forward, B = Backward) while N gives the mean 
number of visits to that class. Here, again EC is Earth-crossing, 
MC is Mars-crossing, SP is short-period and E is Encke-type. 

and P/1998 U3 Jagcr, discussed in Horner et al. 2003), the 
vast majority of such objects are under the control of Jupiter 
prior to a period of cometary behaviour. This is evident in 
Table 6. In both the forward and backward integrations, the 
only classes from which the likelihood of transfer to short- 
period orbits is greater than 0.001 are the other cometary 
classes (E, I, L), the classes in which the object is controlled 
by Jupiter (J, JS, etc.), and the S class. 

Finally, we can use the simulation data to estimate the 
impact rate on the giant planets from Centaurs. In the for- 
ward integration, we find that 144 objects hit Jupiter, 53 
hit Saturn, 5 hit Uranus and a further 5 hit Neptune. In 
the backward integration, these numbers are 135, 48, 5 and 
1 respectively. Given that the estimated population of the 
Centaur region is ~ 44300, then wc expect one impact per 
10 kyr on Jupiter, one per 28 kyr on Saturn, and one per 300 
kyr on Uranus and Neptune. These numbers are likely to be 
underestimates of the impact rate on the planets, given that 
the errors in integration are at their largest when the clone 
is closest to a massive body. 

Of course, all the numbers derived in this Section are de- 
pendent on the assumed flux of one new short-period comet 
every 200 years. If the true flux is higher or lower, then 
the total populations of Centaurs, Earth-crossers and Mars- 
crossers would need to be correspondingly adjusted upwards 
or downwards. Although the flux of new short-period comets 
is perhaps uncertain by a factor of two, it is not uncertain by 
a factor of ten, and so our population estimates are surely 
correct to an order of magnitude. 



6 DYNAMICS AND PHOTOMETRIC 

COLOURS 

Since the Centaurs range over a large area of the outer Solar 
system, it is possible that there could be variations in their 
observable characteristics (such as colour and light-curves) 
as a function of their position. Observations have been car- 
ried out by a number of groups (e.g. Hainaut & Delsanti 
2002; Peixinho et al. 2001; Weintraub et al. 1997; Bauer 
et al. 2003). Although the number of Centaurs with well- 
determined colours is restricted, the situation is rapidly im- 
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Object 


Class 


D 


Ti 


B - V 


Colour 
V - R 


R- I 


2000 EC98 


JU 


F 
B 


0.61 
0.63 


0.854 ± 0.081 


0.47 ± 0.06 
0.466 ± 0.05 


0.50 ± 0.06 
0.439 ± 0.076 


1998 SG35 


JS 


F 
B 


0.67 
0.65 


0.725 ± 0.089 


0.42 ± 0.08 
0.456 ± 0.050 


0.46 ± 0.04 
0.546 ± 0.063 


1999 UG5 


SU 


F 
B 


0.74 
0.85 


0.88 ±0.18 
0.964 ± 0.085 


0.60 ± 0.08 
0.607 ± 0.060 
0.68 ±0.02 


0.72 ±0.13 
0.625 ± 0.042 
0.58 ± 0.02 


Asbolus 


SN 


F 
B 


0.86 
0.75 


0.75 ± 0.04 
0.750 ± 0.040 


0.47 ± 0.04 

0.513 ±0.068 
0.53 ± 0.02 


0.523 ± 0.045 
0.47 ± 0.02 


2001 PT13 


SU 


F 
B 


0.94 
0.87 




0.44 ± 0.03 


0.48 ± 0.02 


2001 BL41 


SU 


F 
B 


0.95 
0.95 




0.53 ± 0.05 


0.56 ± 0.04 


Chiron 


SU 


F 
B 


1.03 
1.07 


0.67 ± 0.06 
0.679 ± 0.039 


0.34 ± 0.03 
0.359 ± 0.027 
0.36 ± 0.10 


0.356 ± 0.037 
0.32 ±0.10 


1999 XX143 


SN 


F 
B 


1.06 
1.38 




0.67 ± 0.07 


0.70 ± 0.06 


Pholus 


SN 


F 
B 


1.28 
1.39 


1.35 
1.19 ±0.10 
1.299 ± 0.099 


0.71 
0.78 ± 0.04 
0.794 ± 0.032 


0.81 
0.814 ±0.056 


1994 TA 


SU 


F 
B 


1.78 
1.52 


1.261 ± 0.139 


0.672 ± 0.080 


0.740 ±0.210 


2002 G09 


UN 


F 
B 


2.93 
3.67 




0.74 ± 0.06 


0.66 ± 0.05 


2000 QC243 


U 


F 
B 


3.18 
3.44 


0.724 ± 0.062 


0.38 ± 0.06 
0.448 ± 0.044 


0.41 ± 0.06 
0.397 ±0.069 


1998 QM107 


UN 


F 
B 


4.87 
5.65 


0.771 ±0.100 
0.730 ± 0.060 


0.474 ± 0.095 
0.520 ± 0.030 
0.63 ±0.12 


0.368 ±0.102 
0.64 ±0.10 


Nessus 


SE 


F 
B 


4.91 
6.40 


0.88 ± 0.07 
1.090 ± 0.040 


0.77 ± 0.05 
0.793 ± 0.041 
0.74 ± 0.08 


0.695 ± 0.066 
0.64 ± 0.07 


Hylonome 


UN 


F 
B 


6.37 
7.30 


0.643 ± 0.082 


0.41 ±0.10 
0.464 ± 0.059 
0.50 ± 0.07 


0.46 ±0.18 
0.490 ±0.122 
0.52 ± 0.06 


2002 GBIO 


UE 


F 

B 


11.1 
13.1 




0.58 ± 0.07 


0.62 ± 0.07 


Chariklo 


U 


F 
B 


10.3 
9.38 


0.802 ± 0.049 


0.47 
0.479 ± 0.029 
0.49 ± 0.02 


0.55 
0.542 ± 0.030 
0.51 ±0.02 


1998 TF35 


UE 


F 
B 


11.5 
10.8 


1.085 ±0.111 


0.65 ± 0.08 
0.697 ± 0.064 


0.66 ± 0.07 
0.651 ±0.119 


2002 GBIO 


UE 


F 
B 


11.1 
13.1 




0.71 ±0.02 


0.65 ± 0.02 



Table 9. List of the Centaurs with colour information. D is the direction of integration and Ti is the half-life in Myr. The colours are 
taken from Hainaut & Delsanti (2002), Boehnhardt et al. (2002) and Bauer et al. (2003). 
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proving. Colours for 19 Centaurs arc given in Table 9. They 
range from some of the reddest objects in the Solar system 
(e.g., Pholus) to much bluer objects (e.g., Chiron). There 
appears to be no correlation between colour and half-life 
within the data set, or indeed between colour and classifi- 
cation (for example the reddest and the bluest objects on 
the list, Pholus and Chiron, are both controlled by the same 
planet at perihelion). 

We expect objects closest to the Sun to have undergone 
some outgassing, and this would lead to a re-surfacing of the 
object, covering the older, darker and redder material with 
fresh material from the interior of the object. Conversely, 
those objects which have not displayed comctary behaviour 
since the formation of the Solar system arc expected to be 
darker and redder. Other factors that could alter the colours 
of the Centaurs include impacts between objects. Collisions 
might activate areas of the surface of the objects, exposing 
fresh material from the interior, manifesting itself as bluer 
colours. It follows that photometric observations of the Cen- 
taurs do have the potential to provide important clues as to 
the dynamical history of the object which cannot be deter- 
mined by integration alone. 

A good example of this is the case of Chiron and Pholus. 
These two objects have similar half-lives, and are controlled 
by the same planet. We can easily calculate that ~ 60% of 
the clones of Chiron become short-period comets in the for- 
ward integrations, against ~ 30% of the clones of Pholus. 
Yet, it is impossible to determine whether either object has 
been a short-period comet in the past, purely by means of 
the integrations. All that we can calculate are percentage 
probabilities. However, it is possible that the fact that Chi- 
ron is the bluest of the Centaurs listed in Table 9 is caused by 
its well-known comctary activity, which may indicate that 
at some time in the past Chiron was a short-period comet. 
Pholus' red colour may similarly indicate it has not yet been 
transferred into an orbit which would activate the surface 
sufficiently for resurfacing to take place. Hence, we may be 
observing the surface of an object which has not been active 
since the birth of the Solar system. 

Bauer et al. (2003) find that the the Centaurs dis- 
play a wider colour distribution than both Jupiter-family 
short-period comets and Edgewortli-Kuiper belt objects. 
The members of the former population arc all active, while 
the members of the latter are not. An intermediate popu- 
lation such as the Centaurs should at least show properties 
varying between one extreme and the other, and this seems 
borne out by the data in Table 9. 



7 CONCLUSIONS 

Detailed simulations of the evolution of the Centaur popu- 
lation under the gravitational influence of the Sun and the 

4 giant planets (Jupiter, Saturn, Uranus and Neptune) have 
been carried out. 23 328 test particles were created by taking 
the orbits of 32 well-known Centaurs and producing 9 clones 
in each of semi-major axis (a), eccentricity (e) and inclina- 
tion (i), giving a total of 729 clones of each object. The clones 
were then integrated in both the forward and backward di- 
rections for a period of 3 Myr. Any clone which reached a 
heliocentric distance of 1000 au was deemed ejected from 
the Solar system and removed from the integration. 



The half-lives for ejection of the Centaurs were calcu- 
lated from the simulation data. These ranged from ~ 540 
kyr for 1996 AR20 (in the forward direction) to ~ 32 Myr 
2000 FZ53 (in the backward direction). The half-life of the 
ensemble of Centaurs was ~ 2.7 Myr, irrespective of the di- 
rection of integration. To analyze the simulation data, we 
exploited a classification scheme introduced by Horner et 
al. (2003), which breaks down the orbits of objects in the 
Solar system into difi^erent classes according to the planet's 
controlling the perihelion and aphelion. The use of the new 
classification scheme allowed us to determine the dynamical 
pathways through which the Centaur population evolves. 
Transfer probabilities were calculated between the different 
classes. These were found to be remarkably independent of 
the object under study, offering the prospect of determining 
possible future histories for objects from the knowledge of 
their current classification. 

Both observations and simulations suggest that one new 
short-period comet is produced every 200 yrs. This is used 
to normalise our simulation results. In this case, the total 
population of the Centaur region (from g > 4 au to Q < 60 
au) with nuclei large enough to provide visible comets (d > 1 
km) is estimated to be ~ 44 300. A flux of 1 new object 
into the Centaur region from the Edgeworth-Kuiper belt 
every 125 yrs is required to maintain the population in a 
steady-state. Additionally, one fresh Earth-crossing object 
is expected to arise from the Centaur region every ~ 880 
yrs. This is both of interest and concern, as large Centaurs 
entering the inner Solar system are likely to fragment with 
production of much dangerous dust and debris. 

In a companion paper (Horner ct al. 2004), individual 
Centaur clones are discussed, showing examples of stable 
resonant behaviour, Kozai instabilities, capture into satellite 
orbits, and evolution into objects spending long time periods 
in the inner Solar system, amongst other things. This gives 
a feel for the rich variety of behaviour of the diverse objects 
known as "Centaurs". 
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